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ABSTRACT 

We present results of the prompt, early, and afterglow optical observations of five 7- 
ray bursts, GRBs 100901A, 100902A, 100905A, 100906A, and 101020A, made with 
the Mobile Astronomical System of TElescope- Robots in Russia (MASTER-II net), 
the 1.5-m telescope of Sierra-Nevada Observatory, and the 2.56-m Nordic Optical 
Telescope. For two sources, GRB 100901A and GRB 100906A, we detected optical 
counterparts and obtained light curves starting before cessation of 7-ray emission, 
at 113 s and 48 s after the trigger, respectively. Observations of GRB 100906A were 
conducted in two polarizing filters. Observations of the other three bursts gave the 
upper limits on the optical fiux; their properties are briefly discussed. More detailed 
analysis of GRB 100901A and GRB 100906A supplemented by Swift data provides the 
following results and indicates a different origin of the prompt optical radiation in the 
two bursts. The light curves patterns and spectral distributions suggest a common 
production site of the prompt optical and high-energy emission in GRB 100901A. 
Results of spectral fits for GRB 100901A in the range from the optical to X-rays 
favor power-law energy distributions and a consistent value of the optical extinction 
in the host galaxy. GRB 100906A produced a smoothly peaking optical light curve 
suggesting that the prompt optical radiation in this GRB originated in a front shock. 
This is supported by a spectral analysis. We have found that the Amati and Ghirlanda 
relations are satisfied for GRB 100906A. An upper limit on the value of the optical 
extinction on the host of GRB 100906A is obtained. 

Key words: telescopes - gamma-ray burst: individual: GRB 100901 A, GRB 100902A, 
GRB 100905A, GRB 100906A, GRB 101020A - gamma-ray burst: general. 
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1 INTRODUCTION 

Since 1997, when the optical radiation of 7-ray bursts 
(GRBs) was first detected, it is kn own that they are the most 
energetic events in the Universe (jKulkarni et al.l[l99 8^). Op- 
tical emission observed hours after a GRB is attributed to 
the so-called afterglow, which is the result of a shock propa- 
gatin g outward in the surrounding media (|Meszaros fc Reed 
1 19971 ). Characteristics of such emission are defined mainly by 
conditions in the interstellar media and the amount of the 
released energy but depend weakly on details of the central 
burst. 

The physics of GRBs and emission mechanisms are not 
completely understood. A further progress calls for more 
observational data and model analysis. An acknowledged 
model is that a GRB is a manifestation of a formation of 
a rotating black hole (or another compact relativistic ob- 
ject) in a course of the gravitational collapse. An engine 
works that converts energy of the collapse into emissions 
of different types, among which a high-energy 7-emission is 
produced for up to several tens of seconds. The afterglow 
emission is detected long afterwards, while the engine is still 
working or not, which is uncertain. In order to understand 
better the details of the process, it is necessary to observe 
the main event itself, at different wavebands, while the en- 
gine is at its most active stage. 

However, to observe prompt optical emission is a chal- 
lenge as a GRB usually lasts no more than several tens 
of seconds. Beginning wit h the promp t optic al observa- 
tions of GRBs in 1998 of lAkerlof et all (|2000l '). successful 
prompt optical detections remain rare. Evidently, two ap- 
proaches can be conceived: to observe extensive sky fields 
waiting for a GRB to occur or to use special robotic tele- 
scopes ready to point anywhere by an alert from an orbital 
7-ray observatory-'alert observations' technique. The Mo- 
bile Astronomical System of TElesc ope- Robots in Russi a 
(MASTER-lfl) use both techniques (|Lipunov et al.ll2010h . 
The present paper is dedicated to the alert MASTER obser- 
vations of five GRBs in Siberia, Ural, and North Caucasus. 

Alert ground-based optical observations of GRBs are 
a new global physical experiment available since the last 
decade. It is made possible thanks to the implementation of 
the global Internet network, powerful personal computers, 
and fast optical CCD-receivers. The challenge is to accom- 
plish quickly the four following steps as early as possible: 

(i) a GRB is d etected by a 7-ray t elescope on board of a 
space craft (Swift fCehrels ct al.'2004'), Fermi l|Atwood et al.l 
I2OO9I ). INTEGRAL ( Winkler ct al. 200i), etc.) 

(ii) after on-board processing is finished, location of a 
GRB is sent to the Gamma-ray bursts Coordinates Network 
(GCN) at NASA. The first two steps take from 10 to 40 
seconds. 

(iii) the burst position is then disseminated to ground- 
based robotic telescopes through the Internet network - in 
about 0.5 s. 

(iv) robotic telescopes are scheduled and pointed to the 
received positions, which takes from 7 to 40 s for moderate- 
size instruments (less than 0.5 m) and from several minutes 



1 MASTER web site: |http://observ.perepIet.ru| 



to hours for 2m- and larger instruments. Whereupon, imag- 
ing is made in optical and IR bands. 

The first Russian robotic telescope MASTER (Mobile 
Astronomical System of Telescope-Robots) came into op- 
eration in 2002 near Moscow thanks to the private fund- 
ing of Moscow association 'Optics0. Const ruction of the 
all-R ussia network MASTER began in 2008 (|Lipunov et al.1 
l2010f) . Presently, telescopes of the MASTER- net are located 
in observatories of Moscow State University (in Kislovodsk), 
Ural State University (in Kourovka), and Irkutsk State Uni- 
versity (in Tunka near Baikal Lake) and the Blagoveschensk 
Pedagogical University (in Blagoveschensk region). These 
observatories span six time zones. Description of the MA S- 
TER 11 telescopes can be found in lKornilov et all (|201ll ). 

In September and October 2010, the MASTER tele- 
scopes made five target observations of GRBs that triggered 
the Swift observatory: GRBs 100901A, 100902A, 100905A, 
100906A, and 101020A. In the present paper, we report these 
observations and make an analysis for two long GRBs, whose 
counterparts were successfully detected. For one of them, 
GRB 100906A, we also use the data obtained with the 1.5-m 
telescope at Sierra-Nevada Observatory (OSN) and 2.56-m 
Nordic Optical Telescope (NOT). 

We begin in Section [2] with a description of the ob- 
servations. Reduction and analysis of the data obtained by 
MASTER are described in Section [S] Preparation of Swift 
light curves and spectra from the data available online is 
also described there. In Section U we turn to the results for 
GRBsl00901A and 100906A: the spectra before T90, the 
late-time spectra and the estimates of the optical extinc- 
tion in the host galaxies, the search of the jet break time 
on the GRB 100906A light curve after 10 000 s. In Section[5] 
we consider the arguments for different sites of the prompt 
optical emission in GRBs 100901A and 100906A. We also 
discuss their spectral evolution, and the spectral character- 
istics of other GRBs, for which optical counterparts were 
not detected. The applicability of the Amati and Ghirlanda 
relations to GRBs 100901A and 100906A is checked. We also 
mention the model of a two-step collapse for the long engine 
activity that can cause a bump on the X-ray and optical 
light curve of GRB 100901A. We summarize our results in 
Section [6l 

Throughout this work, we adopt a standard CDM cos- 
mology with Ho — 70kms~^ Mpc~^, fim ~ 0.27, and 
SIa ~ 0.73. We define the flux density power-law temporal 
and spectral decay indices as Fi, oc t~°'iy~^ , and the photon 
index equals /3 + 1. All errors in the paper are la uncertain- 
ties unless otherwise noted. All Swift data used to reproduce 
light curves and spectra were taken from the Swift Catalog, 
the UK Swift data archival and the Burst Analyser websitelf] 

2 OBSERVATIONS 

Presently the MASTER II net consists of fou r identical, fast 
telescopes located at different sites in Russia l|Lipunov et al.l 

^ http://www.ochkarik.ru/master/ 

http:/ /heasarc. nasa.gov/docs/swift/archive/grb_tabIe/ 
* http:/ /www. swift. ac.uk/swift_portal/ 
^ http:/ /www. swift. ac.uk/burst_analyser/ 



Optical observations of five GRBs 3 



Table 1. Photometry of GRB 100901A by MASTER: R, unfiltered W, I, and V. Times t - TO are the middle times of exposures, which 
are listed in the 2nd, 5th, 8th, and 11th column. Full table can be found online in the electronic issue. 
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4.1381 
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18.08 ±0.45 


0.0315 
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17.00 ±0.09 
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180 
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180 


16.96 ±0.09 
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180 
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180 
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180 


17.61 ± 0.10 
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60 
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180 


16.87 ±0.08 
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180 


17.59 ±0.11 


4.5705 


180 


17.45 ± 0.08 


0.1185 


80 


17.51 ±0.25 


5.6426 


180 


17.04 ±0.09 


5.6451 


180 


17.70 ±0.12 


4.6278 


180 


17.77 ±0.10 
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100 


18.47 ±0.50 
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180 


17.04 ±0.09 
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180 


17.92 ±0.18 
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180 


17.33 ± 0.07 


0.1820 


120 


18.85 ±0.70 


5.7568 


180 


17.05 ±0.09 


7.1718 


180 


17.79 ±0.15 


4.7428 


180 


17.56 ± 0.09 


0.2233 


150 


18.62 ±0.60 
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180 


16.83 ±0.07 








4.8000 


180 


17.45 ± 0.07 


0.2732 


180 


17.94 ±0.30 


5.8706 


180 


16.81 ±0.07 








4.8571 


180 


17.39 ± 0.07 


0.3272 


180 


17.99 ±0.30 


5.9277 


180 


16.92 ±0.07 









l2010l : lKornilov et al.ll20l"lh . Each telescope is a twin-tube in- 
strument with the aperture of 0.4 m and the focal ratio f/2.5. 
Each telescope is equipped with the CCD camera 4 x 4K 
AltaU-16M providing a field of view of 2x2 square degrees. 
A photometric unit provides a set of four filters, which can 
be placed in turn at the optical axis before the CCD camera. 
Depending on the long-term scientific task, the photometric 
units of the twin telescopes are supplied with three filters 
(any of the following: Johnson B, V, R, I, and white glass). 
The fourth position is occupied by polarization filters, ori- 
ented orthogonally in the two tubes. Polarization filters are 
positioned differently in the celestial coordinate system at 
the telescopes of the MASTER net. 

Each observatory operates automatically. Waiting for 
an alert signal, MASTER II monitors the sky searching for 
possible optical transients. Swapping from the monitoring 
to alert observations of GRBs is accomplished by switching 
the filters, focusing, and pointing the telescope to the event 
coordinates received from the GCN. In 20 seconds after an 
alert, the MASTER II observatory is ready to observe a 
GRB. The pointings at the positions of the GRBs considered 
below took from 22 to 47 s. 



2.1 Prompt optical observations of GRB 100901A 

GRB100901A triggered the Swtft Burst Alert Telescope 
(Swift/BAT) at 13:34:10 UT on September 1, 2010. Time 
Tgo, during which 90 per cent of 7-fluence is detected , was 
439 ± 33 s (jlmmler et al.ll2010l : ISakamoto eral1l2010bl ). 

Two MASTER robotic telescopes were pointed to the 
Swift BAT position of GRB 100901A: at Tunka at 103 s 
after the BAT trigger time TO a nd near Blagoveschensk, 
101 s after TO jlvanov et al.ll2010al ). We have got more than 
~ 5 min of optical observations simultaneous to the prompt 
7-emission of the GRB. The source was too low above the 
horizon at Tunka, and it was possible to observe it with only 
one of the optical tubes. The second tube was shadowed by 
the dome, which made impossible measuring polarized light. 
The eastern telescope, near Blagoveschensk, was not oper- 
ational because of the bad weather conditions. An optical 
flare, with maximum at 17.0™ in the MASTER unfiltered 
band, was clearly detected at 426 ± 40 s after the trigger 
time (llvanov et al. 2010bl) . synchronously with the X-ray 
flare dPage fc Immleij|20ld ) and the 7-flare ijSakamoto et al.l 
l2010bl ). 
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Figure 1. GRB 100901A optical light curve obtained by MAS- 
TER (blue points with vertical bars representing errors and hor- 
izontal bars representing exposure times; data are corrected for 
the Galactic extinction) along with the SwiftfBAT 15-350 keV 
flux with 10 s binning (grey dots connected by dot-dashed 
lines) and S'i»ift/XRT 0.3—10 keV unabsorbed flux (black dots; 
iPage fc Immleill2010l'l . Thin-line rectangles show the time inter- 
vals selected for spectral analysis. Color figures can be viewed in 
the electronic version. 



The MASTER telescope at Kourovka was pointed to 
GRB 100901A ~ 5 hours later. Overall, we have been con- 
tinuously observing the GRB during 14 hours at three sites, 
Tunka, Kourovka, and Kislovodsk, in different optical bands: 
unfiltered, R, I, and V. The resulting data are summa- 
rized in Table [T] and presented in Figs. [T] and [51 In Fig. [T] 
Smft/XRT 0.3-10 keV and Swtft/BAT 15-350 keV light 
curves for GRB 100901A are shown as well (for description 
see Section [3. 3p . 

OSN observed the GRB abo ut 11 hours after the trig- 
ger (|Sanchez-Ramirez et al.lboiOl ). The redshift of the host 
galaxy z= 1.408 was obtained at the 8-m Gemini-North tele- 
scope (IChornock et al.ll2010l ). 

2.2 Prompt optical observations of GRB 100902A 

On September the 2nd, the Kourovka MASTER tele- 
scop e was pointed to t he BAT position of GRB 100902A 
( Sak amoto et al.ll20103 ) 104 s after the GRB trigger time 
TO 19:31:54 UT. This GRB flaring activity lasted for 
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Table 3. MASTER upper limits on the optical flux obtained with 
polarizing filter for GRB 100905A. 



t-tr. 



(lirs) 



Figure 2. MASTER optical light curves of GRB 100901A in dif- 
ferent filters: R (3 mag added; red crosses), / (2 mag added; 
magenta dots), V (1 mag subtracted; black squares), and unfil- 
tered (blue dots connected with solid line). Magnitudes are not 
corrected for the Galactic extinction. Horizontal bars show expo- 
sures. 



Table 2. MASTER upper limits on the optical flux obtained 
with the polarizing filter for GRB 100902A. The bottom row is 
coadded from all above. 



^start 
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(UT) 
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19:33:38 


19:33:58 


20 


104.0 


17.0 


19:34:12 


19:34:42 


30 


132.0 


17.3 


19:34:57 


19:35:57 


40 


177.0 


17.3 


19:35:51 


19:36:41 


50 


231.0 


17.4 


19:33:38 


19:36:41 


140 


104.0 


18.0 



about 200 s after the trigger time, and MASTER was 
first to witness its active stage in the optical waveband. 
The imaging was attempted at two polarization angles 
with 20 s exposure times. However, no optical transient 
was fou nd at the XRT position with the upper limit of 
17 mag llKrushinski et aLll2010al). T he final results of our 
photometrv l|Gorbovskov et al.l 201ol ). which are the upper 
limits on the optical flux at two polarization angles, are sum- 
marized in Table [21 The GRB position was observed by the 
OSN at ~ 7 h after the trigger wi th an upper limit on the 
magnitude in /-band: ~ 21.8 mag (jTello et al.ll2010l ). 



2.3 Early optical observations of GRB 100905A 

On September the 5t h, a GRB was de tected by Swift/BAT 
at 15:08:15 UT (Mar shall etld] |2010| ). The Tunka MAS- 
TER telescope was pointed to its position 55 s after 
TO. The imaging was c arried out at two polarization an- 
gles l|lvanov et al.|[2010cl ). No optical transient is found. Up- 
per limits are given in Table |3l The further optical ob- 
servations were carried out by UVOT ( Sie g el & Marshall 
and the 1.5-m telescope of the Observatorio d e Sierra 



Nevada (jSota. de Uearte Postigo fc Castro- Tiradol l2010l'l 
and yie lded no optical counterpart either. iBarthelmv et al] 
l|2010ah obtained that Tqo = 3.4 ± 0.5 s implying that our 
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30 
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40 
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observations can be qualified as 'early', not prompt. This 7- 
ray burst might belong to the class of intermediate duration 
GRBs, separated from the class of long GRBs on statistica l 
grounds by Horvath (1998, 2002): iMukheriee et al.1 ljl998h . 
While underlying physical properties of the subclasses of 
long GRBs appear to be similar (type of progenitor, inter- 
stellar environment), there is an evidence for duration and 
luminosity differences caused by some slight diversity of the 
bursts (jde Ugarte Postigo et al.ii2011i ). 



2.4 Prompt and afterglow optical observations of 
GRB 100906A 

At 13:49:2 7 UT, Septembe r the 6th, GRB 100906A was 
detected (|Markwardt et all I2OIOI ). The BAT /Swift de- 
tected several br i ght pe aks with Tgo — 114.4 ± 1.6 s 
(IBarthelmv et al.l i2010bh . GRB 100906 A has been also 
obse rved by the Korms- Wind starting from 13:49:30 
UT (|Golenetskii et al.ll2010l ). 

Two MASTER telescopes, at Tunka and 
Blagoveschensk, were pointed to the GRB position, 
respectively, 38 s after the BAT trigger time TO and 58 s 
after TO. As a pilot instrument was mounted at the time at 
the Blagoveschensk sit^fl, we present only the data obtained 
by the telescope at Tunka. In Kourovka and Kislovodsk, 
the weather conditions were not suitable for observations. 

A bright optical transient wa s localized by MA STER at 
the 5ot/£/UVOT position (,Markwardt et al.lf201ol ). 13 mag 
at maximum. During the first hour after the trigger time, 
24 images in unfiltered light at two polarization angles were 
produced by the Tunka telescope. Fig. [3] shows the early 
observations of GRB 100906A in optical, 0.3-10 keV, and 
15-150 keV energy range (also Fig. |4]and Table|4]). A movie 
of the optical burst of GRB 100906A as seen at Tunka is 
available on the welQ. 

We calculate the relative difference between the signals 
coming from the two polarization filters, for the time inter- 
val 100 — 10^ s after TO. It is less than the relative mea- 
suring accuracy of filters 2 per cent, estimated through the 
observations of standard stars. Unfortunately, technical lim- 
itations do not allow us to draw a decisive conclusion about 
the polarization degree of GRB 100906A (see § 15.11 and the 
Appendix). 

The 1.5m OSN telescope at Sierra Nevada observed the 
GRB few hours after TO (Fig. [S] and Tabled. Calibration 



^ Since December 2010, the Blagoveschensk telescope is in a full 
operating mode. 

^ http: / /master. sai.msu.ru/static/GRB/grbl00906.avi 
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Table 4. Photometry data for GRB 100906A by MASTER in the 
unfiltered band with two orthogonal polarizing filters. Full table 
can be found online in the electronic issue. 
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Figure 4. MASTER light curve of GRB 100901A 
izations, unfiltered band (blue dots and red crosses) 
corrected for the Galactic extinction. 



in two polar- 
Data are not 




100 1000 

Figure 3. GRB 100906A optical light curve obtained by MAS- 
TER with one polarization filter (blue dots with 20 per cent error 
bars resulted due to uncertainty of magnitude-flux conversion, 
see § 13.21 data are corrected for the Galactic extinction) along 
with the Swift/BAT 15-350 keV flux with 1 s binning (grey dots 
connected with dot-dashed line; data with negative lower limits 
are not shown) and t he 5mft/XRT 0.3-10 keV unab sorbed-fiux 
light curve (black dots: lBeardmore fc Markwardtl2010l ). Thin-line 
rectangles show the time intervals selected for spectral analysis. 



for BVIR bands is done using 49 USNO stars, and for U- 
band using a Landolt star observation. Aperture photome- 
try is carried out using phot as implemented in IRAF with a 
radius equal to the seeing. The error of the afterglow magni- 
tude was obtained adding quadratically the statistical error 
given by phot and the zero point error determined with the 
49 field reference stars. In addition. Table |B] shows the late 
photometry obtained by the 2.56-m Nordic Optical Tele- 
scope (NOT). 

The redshift of the host galaxy z = 1.727 
was obtained with the Ge mini-North tele- 
scope (jTanvir. Wiersema fc LevanlbOld '). 
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Figure 5. OSN and NOT telescopes i?-band light curve of 
GRB 100906A (red bars and green triangles with small error bars, 
respectively). MASTER unfiltered band, one polarization light 
curve is plotted by blue crosses before 5 h. All data are not cor- 
rected for the Galactic extinction. 



Table 5. !7Syi?7-photometry for GRB 100906A obtained by the 
OSN telescope. Reported times are middle times of exposures. 
Dots substitute 50 observations that can be found online in the 
electronic version. 
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Magnitude 


R 


7.920835 


36 


18.53 ±0.12 


R 


7.935742 


36 


18.72 ±0.12 


R 


8.895094 


60 


18.79 ±0.06 


R 


14.091390 


180 


19.49 ± 0.06 


R 


14.217964 


180 


19.50 ±0.06 


R 


14.540464 


180 


19.55 ±0.06 


R 


14.784353 


180 


19.76 ±0.08 


R 


14.928057 


180 


19.86 ±0.08 


R 


38.528342 


2x300 


> 21.5 


V 


14.141108 


180 


19.91 ±0.20 


V 


38.402228 


300 


> 21.0 


B 


14.193888 


300 


20.65 ±0.22 


B 


38.318343 


5x300 


> 21.0 


I 


14.278609 


2x180 


19.46 ±0.12 


U 


14.571669 


2x900 


> 17.5 


U 


38.234160 


300 


> 16.5 
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Table 6. R-magnitudes of GRB 100906A obtained by NOT. 



tstart - TO i<;nd " TO Exp. time R 
(h) (h^ (s) (mag) 

35.031028 36.338862 3 x 600 22.15 ±0.05 
62.828250 63.078251 900 23.10 ±0.09 

84.046943 84.763888 2 x 900 23.56 ± 0.06 



Table 7. Upper limits on the optical flux in V-band for 
GRB 101020A obtained for coadded images. 



tatart — TO tmiddic Exp. time V Coadded 

(s) (s) (s) (mag) number 

106 116 20+20 14.9 1+1 

472 2836 4300 18.5 25 



2.5 Limits on the prompt and afterglow optical 
flux of GRB 101020A 

On October the 10th, at 23:40:41 UT, the gwft/BAT 
trlRKered and located GRB 101020A l|Saxton et al.ll2010l ). 
ISakamoto erall (|2010al ) found Tgo = 175 ± 28 s. MASTER 
telescope at Kourovka was pointed to the GRB 101020A po- 
sition 106 s after TQ at a lar ge zenith distance of 80 de- 
grees (jKrushinski et al.ll2010bl ld). It was observed with two 
optical tubes. No optical transient is detected with the up- 
per limits in V-band reported in Table [71 where magnitudes 
are obtained from coadded images. 



3 DATA REDUCTION AND ANALYSIS FOR 
GRB 100901A AND GRB 100906A 

3.1 MASTER astrometric and photometric 
calibration 

Each MASTER telescope is equipped with a set of John- 
son/Bessell filters and two linear polarizing laminated films. 
For every image, astrometric calibration using imcoords of 
IRAF v2.14 is performed. Aperture ph otometry is done with 
the phot of IRAF package (Tody 1993). Standard stars from 
13.5™ to 17™ are selected fro m the Sloan Digital Sky Survey 
Data Release 7 (SDSS-DR 7) (|Abazaiian et al.ll2009l '). Influ- 
ence of atmospheric transparency variations on instrumental 
stellar m agnitudes is eliminated f ollowing an algorithm de- 
scribed in lEverett &: Howelll (|200ll ). Stars with variations ex- 
ceeding threefold the errors calculated by IRAF from signal- 
to-noise information are expelled from the list, and the re- 
maining standard stars are related to SDSS-DR7 data as 
suggested by Lupton (20050: 

B = 5 + 0.1884 (1^-3) + 0.1313, (1) 
V = g~ 0.5784 {g - r) - 0.0038 , 
R = r - 0.2936 (r-i)- 0.1439, 
I = i~ 0.3780 {i~z)~ 0.3974 , 

where BVRI are the Johnson system and i, r, g, and z 
are the Sloan Survey photometric system (all in the Vega 

^ http:/ /www. sdss.org/dr7/algorithms/sdssUBVRITransform. html 



system). Resulting lists of standard stars are given in Ta- 
ble m for GRB 100901A and in Table |9] for GRB 100906A. 
Associated variances'' between observed and catalogue val- 
ues after applying the above relations in different filters are 
0.005 — 0.008™. In our observations, standard deviations of 
light curves of the standard stars are greater, 0.02 — 0.1 mag 
depending on the filter. For images obtained in the unfil- 
tered band, we use a zero-point correction to the unfiltered 
magnitudes obtained at the Kislovodsk MASTER telescope. 
Flux calibration of the unfiltered magnitudes at Kislovodsk 
is described in Section [3.21 

The GRB magnitude errors are defined by the standard 
deviation of the brightness variations of the standard stars 
of similar magnitude in the case of high signal-to-noise ra- 
tio. At the low signal-to-noise ratio, we lack standard stars 
of such magnitudes. To estimate observational errors, we 
take the theoretical limit for errors following the basic CCD 
equation. This limit is much greater (up to 2 mag at mag 
> 18) than expected amplitude variations of standard stars 
at such signal-to-noise level. 

For each of the standard stars, we evaluate the differ- 
ence between its average magnitudes in both orthogonal po- 
larizations. Assuming that the light of standard stars should 
not be polarized, the standard stars with outlier polariza- 
tion are expelled from the list. Afterwards, corrections are 
applied to match average magnitudes of the standard stars 
in two polarizations. The relative accuracy of observations 
with two filters is approximately 2 per cent, basing on ob- 
servations of standard stars. 

3.2 MASTER flux calibration 

To convert magnitudes BVRI into absolute values, we use 
the zero-magni tude flux densities of th e Landolt photom- 
etry (see, e.g.. Pic kles fc Depagii3l2O10l ). To convert CCD 
photometry into absolute fluxes without color filter inter- 
posed, that is, in a wide waveband (see Fig. [SJ, we use the 
Pogson formula relating stellar magnitudes and CCD band- 
integrated flux with a zero- ma gnitude flux c alculated using 
the Vega spectrum as given by iHavea ((1985t ). 

The standard stars are designated with the W and P 
magnitudes: obtained from observations in the wide wave- 
band and in the wide waveband plus polarizer. Both are 
approximated by the formula 0.8 R + 0.2 B. This formula is 
arrived at by compilation of color diagrams W — X versus 
B — R, where X is a combination of stellar magnitudes, 
X = aR + bB. In a general case, points group around 
a straight oblique line. The zero-slope of the relation is 
achieved by fitting for the parameters a and b. It is found 
that for the white light the above representation of the wide- 
band magnitudes holds with an accuracy of 4 per cent for 
W— magnitudes and 2 per cent, for P— magnitudes. Unfortu- 
nately, for a polarizer that we used in the observations there 
was a long- wavelength interval with unspecified transmission 
efficiency function. Simulating polarizer transmission func- 
tion in different ways and folding it with different spectra 
(Vega's, cx , and oc v^), we have estimated the maximum 
additional uncertainty of ~ 20 per cent. Thus we add 20 per 
cent of uncertainty when we convert P to absolute fluxes. 

Observed magnitude W of an investigated source is con- 
verted into the absolute flux values using the Pogson equa- 
tion 
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(2) 



3.4 Spectral fitting 



where is the calculated Vega flux in the CCD spectral 
band: 



fVcga(A) W{\) dX = 1.33 X IQ-'^ergcm-^ s"' . (3) 



In the same way, magnitudes P are converted to flux values 
using eq. ([2]) with the zero-magnitude flux 



= J Fvega(A)P(A)dA = 1.16 X lO"'^ erg cm"^s"V 



(4) 



In the above expressions, W{X) and P(A) are, respectively, 
the normalized CCD-response function and its normalized 
convolution with the transmission efliciency of the polarizer 
(see the right panel of Fig. [6|. Of course, if one observes 
a source with partly (linearly) polarized light through such 
polarization filter, one needs a full information on the degree 
of polarization and its angle to derive the total flux. For an 
unpolarized source, or for a source with equal signals from 
two polarizers (the case of GRB 100906A), we can safely use 
formula ((4)). For details we refer to Appendix. 

To obtain GRB's flux density at the wavelength of the 
CCD-response maximum, 5500 A, we divide flux -Forb or 
^GRB by an effective frequency interval Aucs of the corre- 
sponding response function. Naturally, the value of AvcB 
depends on a particular spectral distribution. The effec- 
tive frequency interval of the CCD-response function is 
~ 3.9 X 10^* Hz, calculated as the mean of values Aucs for 
the Vega spectrum, white light, and power-law spectrum 
with —1 slope, all falling inside a 15 per cent uncertainty 
interval. The effective frequency interval of the CCD plus 
polarizer, calculated in the same way, is ~ 3.2 x 10^* Hz 
with a 10 per cent uncertainty interval. Bearing in mind 
that a real spectrum might be something different from the 
three considered in such simple analysis, we assume approx- 
imately a 20 per cent accuracy for converting from flux units 
to spectral flux density units. 



3.3 Swift BAT and XRT light curves 

In order to build flux light curves observed by 
Swift/BAT (Figs. [Hand [3|, standard sets of GRB products 
are obtained using the ftool batgrbproduct within HEA- 
SOFT v6.10, and the count light curves in 15-350 keV are 
retrieved. The 15 — 150 keV spectra averaged over the total - 
fluence time interval are analysed in XSPEC (|Arnaudlll99"6l ') 
using the power- law approximation. The derived spectral pa- 
rameters are used to calculate 15 — 350 keV model flux and to 
obtain counts to flux conversion factors in the 15 — 350 keV 
range. Thus we apply constant conversion factors 1.07 x 10~® 
erg cm"^ counts"^ for GRB 100901A and 9.17 x lO"'^ erg 
cm"^ counts"^ for GRB 100906A. 

The Swift/XRT 0.3-10 keV unabsorbed light curves are 
available from the Swift Burst Analyser. These light curves 
are pro duced with inform ation on spectral evolution in- 
volved (|Evans et al.ll2ofol ). The late-time spectra of GRBs, 
available from the Swift Spectra Repository, are the source 
of estimates of the column densities for the soft X-ray ab- 
sorption. 



In order to build simultaneous wide-energy spectra, we select 
time intervals, which are covered by the optical observations. 
The prompt stage of GRB 100901A (before Tgo) contains 
three suitable MASTER exposure intervals: t -TO = 193- 
233, 312-372, and 387-467 s, where TO is the trigger time. 
There are two MASTER time intervals during the prompt 
stage of GRB 100906A overlapping with the observations by 
Swift/XRT, 73.8 - 83.8 s and 105 - 125 s (or very close, in 
the flrst case). 

We download the X-ray time-sliced spectra, con- 
structed for the time intervals of the optical expo sures, 
from the Swift spectra repository (|Evans et al.ll2009l ). For 
GRB 100906A we also use the resu lts obtained by the 
Konus- Wind (jColenetskii et al] I2OI0I '). BAT spectra and 
response matrices are prepared as described in the Data 
Analysis Guideline^ They are extracted for the selected 
time intervals and corrected using the tasks batbinevt, 
batupdatephakw, batphasyserr, and a response matrix gen- 
erated by the batdrmgen task. Then the spectra are fltted in 
XSPEC vl2.6 using models powerlaw (power law), cutoff pi 
(cut-off power law), or grbm (Band function), in combina- 
tion with absorption components. When simultaneous fits 
for the 0.3 — 10 keV and 15 — 150 keV and/or optical data 
are performed, XRT spectra are rebinned using grppha to 
ensure that there are at least 20 counts per bin, which is 
needed for joint minimization of x^- 

To account for the soft X-ray absorption, we use XSPEC 
model component phabs (Galactic absorption) and zphabs 
(absorption in the host galaxy). When we fit the parame- 
ters of the optical extinction in the GRB host, we use the 
XSPEC model zdust, which is the extinction by the dust 
grains. In the fits, whose results are presented in the tables 
in the next section, we assume the Small Magellanic Cloud 



|Pei| 


1992h foUowins 


I2OOI 


1: ISchadv et al.l 



(|2010l ). We discuss variations due to the different extinction 
laws in the text. Correction for the reddening in the Galaxy 
is done beforehand, using the values calculated using the 
NED Galact ic Extincti on calculatoiQ based on a work of 
ISchlegel. Finkbeiner fc Davisi (|l9S?3 ). 

To analyse the late-time spectra, we consider the optical 
observations after IOT90. First, we check that big variations 
of X-ray flux are absent in a time interval: we request that 
the flux density at 10 keV varies less by a factor of 3. The 
X-ray data is rebinned to a minimum of 10 (20, if possible) 
counts per bin (by grppha). The optical data is averaged over 
the time interval and transformed by the task f lx2xsp. We 
check that the number of degrees of freedom in the fit by the 
XSPEC model phabs*zphabs*zdust*powerlaw is more than 
5. Thus selected time intervals are depicted by rectangles in 
Figs. [H [31 and [11] 



^ |http;/ /swift. gsfc.nasa.gov/docs/swift/analysis/threads/batgrbproductthread.htn 
http:/ /ned. ipac.caltech.edu/forms/calculator. html 
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Table 8. Stars from SDSS-DR7 catalogue used as standard stars for GRB 100901A. Columns 3-7 contain 
SDSS ugriz magnitudes (Vega system). Column 8 contains average unfiltercd W-band MASTER magnitude. 
Pull table can be found online in the electronic issue. 



Source 


SDSS 


u 


a 


r 


i 


z 


W 


RA J2000 


DEC J2000 


ID 


name 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(d) 


(d) 


1 


J014906. 68+224024.7 


19.246 


17.05 


15.96 


15.462 


15.161 


15.95 


27.277834 


22.673535 


2 


J014904.39+224045.3 


15.841 


14.609 


14.143 


13.953 


13.897 


14.21 


27.26829 


22.679256 



Table 9. Standard stars used for the photometry of GRB 100906A. Column 3 contains average unfiltered 
P-band MASTER magnitude measured with one of polarization filters. Full table can be found online in the 
electronic issue. 



Source 
ID 


USNO-A2.0 

name 


(mag) 


Standard deviation 
of light curve (mag) 


(Pj) - (P„) 
(mag) 


RA J2000 
(d) 


DEC J2000 
(d) 


(1) 
2 


(GRB 100906A) 
1425-02624804 


15.70 


1.14 
0.04 


0.01 
-0.04 


28.68413 
28.69122 


55.63050 
55.63834 



4 RESULTS 



4.1 Spectra of GRB 100901 A and GRB 100906A 
before Tgo 

For the three early spectra of GRB 100901A (Fig. [7)1, fits 
with absorbed power-law are satisfactory (Fits 100901.1, 
100901.2, and 100901.3 in Table[TD]). Other spectral models, 
cut-off power law and the Band function, do not improve 
the fitting statistics. 

For GRB 100906A, we do not include the host-galaxy 
absorption model component, as its spectral distributions 
from the optical to X-rays (Fig. [S]) cannot apparently be 
fitted by a single power law. Generally, for a variable com- 
plex spectrum of an early GRB emission, it is impossible 
to estimate the extinction in the optical band, as one can 
use no a-priori information about the spectral distribution. 
Table [TT] contains the parameters of the fits by the Band 
model, together with the resulting values of for the alter- 
native models with fewer free parameters, which give worse 
results. We also calculate the slopes of the photon spectrum 
between the MASTER points and the best-fit 3 keV flux 
densities (column 3 in Table [TT|) . We remind that the opti- 
cal points are corrected for the Galactic extinction. 

Apparently, the peak of the spectral energy distribution 
(SED) of GRB 100906A is located somewhere around 0.3- 
30 keV at times 80-120 s and moves gradually to the soft 
energies. The SED-peak energy i?peak for the Band spec- 
trum is (— Qb + 2) X Eo, where qb is the photon index of 
the lower-energy part of the Band function, and Eo is the 
characteristic energy for the Band function (Table [TT|) . We 
obtain that -Epcak shifts from about 30 to 4 keV between the 
two spectral distributions presented in Fig. [S] The higher 
energy part of the spectrum of GRB 100906A (> 30 keV) 
can be fitted by a single power law. This is consistent with 
the Konus- Wind results (Fig. [S] dotted line). According to 
iGolenetskii et all (2010), the spectrum of the second burst- 
ing episode of GRB 100906 A observed by Komis- Wind in 
7-rays from 98.304 to 122.880 s is best fitted in the 20 keV- 
2 MeV range by a power-law model with the photon index 



4.2 Host galaxy extinction of GRB 100901A 

We fit several late-time spectra comprised of the XRT/Suii/i 
data and simultaneous optical points available from MAS- 
TER. The optical points are the time-averaged unfiltered W 
magnitudes as follows: 18.05±0.03, 17.78±0.02, 17.61±0.04, 
and 19.30 ± 0.02 mag, corrected for the Galactic extinctio n 
Av = 0.327 (NED; ISchlegel. Finkbeiner fc DavisllTggi '). 
successively for the fits in Table 1121 As well as for the 
prompt spectra, we can approximate the late-time spectra 
of GRB 100901A by an absorbed power-law model. 

Averaging results of the late-time fits (Table I12|) . we 
obtain A^^ = 0.4 ± 0.2 and iV^f' = (2.2 ± 0.7) x 10^^ cm'^ 
Averaging results of the three early fits ( Table [T0||. we obtain 
yl'^* = 0.8 ±0.1 and iV^f* = (5.1 ±0.6) x lO^^cm"^ There is 
an apparent change in values of A^^* and Ay*" between the 
prompt-time and late-time fit parameters. It is likely that 
the prompt-time value of A*'^* is biased because of the spec- 
tral evolution and, probably, because the prompt spectral 
shape is not pr ecisel y powe r-law in the r ange from the op- 
tical t o X-rays. iButler fc Kocevskil (|2007| ) and lEvans et al.l 
(|2010h argue that absorption derived from early-time XRT 
data, when a strong spectral evolution is present, can be 
misleading. 

The late-time spectrum compiled at the yiRT/ Swift 
spectrum repository (UK Swift Science Data Centre) for 
4000-33229 s yields the value of the intrinsic A^^f* = (2.9 ± 
0.9) X lO'^^cm"^ (90 per cent uncertainty interval), which is 
in a very good agreement with our fit results for the late- 
time sliced spectra. This indicates that a power law approx- 
imation for the X-ray spectrum of GRB 100901a is indeed a 
good approximation for any late time interval. 

We have also tried other types of extinction laws in 
the GRB host. In addition to the fits, presented in Ta- 
ble 1121 we performed fits with different extinction laws in 
the host galaxy: that of the LMC (Large Magellanic Cloud) 
and MW (Milky Way). For them we applied different values 
of Rv (jPeiii 19921 , see below) . The quality of the fits remains 
the same and the values of A^h as well, because the choice of 
the optical extinction does not affect strongly the amount of 
the soft X-rays absorption. We have found that depending 
on the extinction law applied, for three late-time intervals, 
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Wavelength (A) Wavelength (S) 

Figure 6. In the left panel the MASTER filter transmission efficiencies are plotted along with those of the CCD and polarizer. In the 
right panel, functions W{X) and P(X) are shown, which are the bandpasses of the CCD and of the CCD plus polarizer, respectively. 
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Figure 7. Spectrum of GRB 100901A for two time intervals at t < Tqq. Op tical flux density obtained by MASTER, corrected for 
the Galactic extinction Ay = 0.327 (NED; ISchlegel, Finkbeiner &: Davislll998h . is shown by the single left point, whose horizontal bar 
corresponds to the MASTER unfiltered effective frequency interval. Spectra in 0.3—10 and 15—150 keV are made with the Swift BAT 
and XRT data. Best-fitting absorbed power laws are shown by the red lines. Their spectral parameters are described in Table [TO] as Fit 
100901.2 for 312 - 372 s and Fit 100901.3, for 387 - 467 s). 




E (keV) E (keV) 

Figure 8. Spectr um of GRB 100906A for two time intervals at t < Tqq. Optical points are corrected for the Galactic extinction 
Ay = 1.194 (NED; ISchlegel. Finkbeiner fc Davis|[l998l '). In the left panel, we use the MASTER observation at 73.8-83.8 s. Best-fitting 
absorbed Band functions are shown by the red lines (Fit 100906.1 and 100906.2 in Table [TOt . No attempt has been made to estimate 
the optical extinction in the GRB host galaxy. The brown dashed line depicts the unabsorbed low-energy part, and the orange dot- 
dashed line, the unabsorbed high-energy part. Dotted power law represents observations of Korms- Wind from 98.304 to 122.880 s in 
20 keV-2 MeV with a correct slope and a roughly estimated flux. 
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Table 10. Spectral parameters for GRB 100901A at the prompt stage resulted from the simultaneous fitting of the absorbed power law 
in XSPEC (phabs*zphab s*zdust*po«erlaw) to the MASTER, XRT, and BAT data (Fig.O. Galactic absorption column density is fixed 
at Afjj = 7.1 X 102''cm~2 (jKalberla et al."2005'). Assumed host galaxy extinction curve is SMC type and the total-to-selective extinction is 
Ry = 2.93. The 4th and 5th columns give the fitted intrinsic absorption column density and optical extinction at z = 1.408, respectively. 



Fit ID 


t - TO (s) 


Photon index 


Af^"t (10^2 cm-2) 


/tint 


xVdof 


100901.1 
100901.2 
100901.3 


193 - 233 
312 - 372 
387 - 467 


1.60 ±0.08 
1.53 ±0.05 
1.49 ± 0.02 


0.6 ±0.2 
0.4 ±0.1 
0.55 ±0.07 


1 1+0.6 

0.951°;^ 
0.75±0.1 


73.4/94 
138/134 
216/219 



Table 11. Spectral parameters for GRB 100906A at the prompt stage resulted from the fitting of the absorbed Band func- 
tion in XSPEC (phabs*2phabs*grbm) to the XRT and BAT data (Fig. [S]!. Galactic absorption column density is fixed at A^h = 
2.2 X 10 cm" -2 jKalberla et al.ll2005h . Column 3 gives the photon index between the optical and modeled 3 keV flux density; next 
columns: the photon indices and characteristic energy of the Band function, intrinsic absorption at z = 1.727. In the last three columns 
statistics of the fits shown in Fig. |8] are typed in bold, statistics of fits with the alternative models are typed in normal font. 



Fit 
ID 


t - TO 


MASTER-XRT 
2.2 eV-3 keV 


Photon indices 
of Band function 


Eo 
(keV) 


(1022 cm~2) 


xVdof 
Band 


xVdof 
power law 


xVdof 
cut-off pow 


100906.1 
100906.2 


84 - 94'l' 
105 - 125 


1.50 ±0.06 
1.24 ±0.06 


1.2 ±0.2 2.18+o;27 
-0.22+°;2* 2.70 ± 0.03 


9 q+0.4 


1+°-^ 

n+0.23 
'^-o.o 


112/110 
149/139 


141/112 
314/141 


121/111 
254/140 



t MASTER observation was carried out in an slightly earlier time interval, 73.8-83.8 s. 



the average total extinction = Rv E{B — V)'"* varies 
as: (A^^) = 0.4 ± 0.2 (SMC, Rv = 2.93), {A^') = 0.5 ± 0.2 
(LMC, Rv = 3.16), and = 0.7±0.2 (MW, Rv = 3.08). 

Basing on the 'late-time' values of the optical extinction 
in the GRB host we calculate the ratio N^*" /Ay^. It equals 
(5 ± 2) X 102^cm"2 (for the SMC type of the extinction 
in the GRB host), (4 ± 2) x 102icm-2 (LMC) and (3 ± 
1) X 102^cm"2 (MW). These are not far from the ratios 
empirically found for the Milky Way and Magellanic Clouds: 
~ 4 X 102i ( SMC), ~ 3.5 X 102i (LMC), and ~ 2 x 102icm-2 
(MW) (e.g.. lSchadv et al.ll2010l '). Thus, independently of the 
extinction law used, the resulting metal-to-dust ratio in the 
host of GRB 100901A is comparable to the level measured 
for the Milky Way or Magellanic Clouds. 



4.3 Host galaxy extinction of GRB 100906 A 

Three late-time spectral distributions are shown in Figs. [9] 
and llOl The parameters of the fits are listed in Table [T3l The 
optical points are the time-averaged unfiltered P magnitudes 
by MASTER: 17.54 ± 0.05 for 6000 - 7000 s, 18.7 ± 0.3 for 
10500-11000 s; OSN R = 18.70±0.02 for for 28000-30400 s. 
The optical point in the 7?-band at 50000-54000 s is ob- 
tained by averaging the OSN data (i? = 19.7 ± 0.3) and 
V, B, I points can be found in Table O Before the fit- 
ting, the optical data are corrected for the Galactic ex- 
tinction calculated in the NED Extragalactic Calculator: 
Av = 1.194, Ab = 1.554, Ar = 0.963, and Ai = 0.699 
(NED; ISchlegel. Finkbeiner fc Davisll 19981 ). 

One can see that the optical data at ~ 14 h (Fig. llOp 
is best fitted by an absorbed broken power-law with a peak 
in the optical range. However, we cannot consider this peak 
robust because the optical points are non-simultaneous and 
there is a strong variation observed in the i?-band around 
14 h (see Fig.[5|. 

The parameters of the Fit 100906.3 (the first three lines 
in Table [T3)) are obtained by simultaneously fitting the data 



for the three time intervals with the spectral slope and nor- 
malization left independent and tied E{B — V)'"' (that is, 
Av) and A^^*. This was done in an attempt to increase 
statistics to obtain positive Icr interval for Ay^. Neverthe- 
less, we can report only the 2a upper limit. Setting other 
types of the extinction law in the host galaxy does not 
change its value. 

The spectrum of GRB 100906A at 10-50 ks inter- 
val from the XRT spectrum repository (at the UK 
Swift Science Data Centre) yields the value of intrinsic 
N^^ = (8.5 ± 3) X 102i cm"2 (with 90 per cent uncertainty 
interval). The values of TV^f* in Table [TS] with la errors are 
consistent with this value. If we take the lower 2a limit of the 
XRT repository value, Af|f* = 4.8 x 102^ cm"2, and the up- 
per 2a limit on Ay^ obtained for different extinction laws, we 
arrive at the lower 2a limit on Nl^^/A'y^ « 1.1 x 1022. This 
limit is at least twice the typical ratios, empirically defined 
for the Milky Way and Magellanic Cl ouds (see This 
ratio is consistent with the results of ISchadv et al.l (|2010l ) 
who studied 28 GRBs and found values of Afg'V^v* typi- 
cally higher than those for the Milky Way and Magellanic 
Clouds. 

We conclude that for GRB 100906A our fitting provides 
only the upper limit on Ay^, and that there are indications 
that this GRB has high N^^^/Ay^ ratio, which can be due 
to a dust destruction or the intrinsics properties of the host 
galaxy. 



4.4 Jet break time of GRB 100906A 

We construct broken power-law fits to the light curves of 
GRB 100906A after 10 ks. Table fTD presents the parameters 
of the fits. A break on the X-ray light curve is found around 
~ 50 ks. The parameters of the 0.3 — 10 keV light curve 
can be reconciled with the "closure relations" between the 
temporal and spectral indices {a and /3 in F,^ oc t~"u~^) 
in the framework of the fireball model in the following way 
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Table 12. Spectral parameters for GRB 100901A at late time intervals. XSPEC fitting of the absorbed power law to the MASTER and 
XRT/ Swift data is performed. Frozen fit parameters are the same as for the fits in Table [Tol 



Fit ID t-TO(s) Photon index N^^ {10^^ cm'^) A'^* X^/dof 



100901.4 


5000 - 


6500 


2.1±0.1 


0.3 ±0.1 




15.8/14 


100901.5 


15000 - 


17500 


2.0 ±0.1 


i+°-2 


o.2i°:i 


38.6/38 


100901.6 


27000 - 


28500 


2.1 ±0.1 


0.3 ±0.2 




36.6/35 


100901.7 


102800 - 


121500 


2 lo+° °i 


0.28 ±0.25 




13.0/10 



Table 13. Spectral parameters for GRB 100906A at several late time intervals (Figs. [9land IIOII . XSPEC fitting of the absorbed power 
law or broken po wer law to the MAS TER, OSN, and XRT/Swift data is performed. Galactic absorption column density is Nn = 
2.21 X lO^^cm"^ llKalberla et alj[2005l 'l. The host galaxy extinction curve is that of SMC (Small Magellanic Cloud) and the total-to- 
selective extinction is taken Ry = 2.93. The columns designation is the same as for Table [121 with an addition of the lower photon index 
(left) and the break energy in eV. 



Fit ID 


t - TO (s) Photon index Bbrcak (eV) 


Photon index 


TVjj^t (10^2 cm-2) 


Aint 


xVdof 


100906.3 
100906.4 


6000 - 7000 
10500 - 11000 
28000 - 30400 

45000 -65000 -1.3 ±0.7 2.5 ± 0.3 


1.96 ±0.05 

1.9 ±0.1 
2.01 ±0.02 
2.34 ±0.04 


25+°-55 


< 0.4 

0.35 ± 0.09 


36.6/42 
9.78/10 



(see table 1 inlRacus n et aP (l2009f). references therein and 
revie ws of iMeszaroj |2002| : IZhang fc Meszaroj |2004| : iPiranI 
[iooi). Before the break, for the ISM environment and the 
'slow cooling' regime and using /3xrt = 1.0 ± 0.3 we get 
Q = 3/3/2 = 1.5 ± 0.45 (c.f. a = 1.68 ± 0.02 in Table [Til 
the 1st row, before the break). The electron spectral index is 
p = 2 /3 + 1 = 3.0 ± 0.6. After the break, applying a uniform 
spreading jet scenario,we arrive aia — 2(5 + 1 — 3.2 ± 0.4. 
The electron spectral index has the same functional form 
P = {p — l)/2. The value of /3 is almost the same, (/3xrt = 
1.1 ± 0.2), and we arrive at p = 3.2 ± 0.4. Thus the closure 
relations can be satisfied for a consistent value of p. The 
relations applied imply that Vm < < Vc, where Vm and 
are the characteristic synchrotron and cooling frequencies, 
respectively. 

We have tried to find out whether this break is achro- 
matic. The OSN data R show that around 52 ks there is an 
abruptness of the R light curve (the fit parameters in the 
row 2 of Table [H]). But the time sampling of OSN R data 
alone is not good enough to clear up the behavior around 
this time. If we combine the OSN and NOT data (the 3rd 
row in Table I14p , we apparently do not get an achromatic 
peak as the tentative break on the R light curve shifts to 
about 35 ks. 

The slope of the R light curve (OSN+NOT) before the 
break is obtained by the data, crowded around 30 ks, and 
might be misleading. If we suggest that the slope of the R 
light curve before the break equals to the slope of the MAS- 
TER light curve for the earlier data (blue crosses in Fig . llip . 
that is 1.08 (also reported by iKuvshinov et aL I I2OI0I ). the 
break time shifts to 39.0l°'^ ks (the temporal slope after 
the break remains the same ^ 2.2 as in the 4th row of Ta- 
ble [H). 

We conclude that the X-ray data alone indicates a jet 
break that can satisfy the closure relations. The optical data 
do not yield confidence about the achromatic character of 
the break. For the discussion of Amati and Ghirlanda rela- 
tions, see i) 15.41 
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Figure 9. Spectral distributions of GRB 100906A at three time 
intervals: 6000-7000, 10500-11000 and 28000-30400 s (points, 
squares and crosses with bars, respectively). The optical data are: 
the MASTER P = 17.54±0.05 for 6000-7000 s, P = 18.7±0.3 for 
10500-11000 s and OSN R = 18.70±0.02 for for 28000-30400 s, 
additionally corrected for the Galactic extinction Ay = 1. 194 and 
Afi = 0.963 (NED; ISchlegel. Finkbeiner fc David Hggj) . Lines 
show the best-fitting absorbed power laws, whose parameters are 
listed as Fit 100906.3 in Table [Tsl 



5 DISCUSSION 

5.1 Optical emission sites of GRBs 100901A and 
100906A 

The prompt behavior of the two closely studied GRBs is 
quite distinct. For GRB 100901A, we see a similarity of light 
curves in the optical, 7, and X-rays at around 300-500 s after 
the trigger. For GRB 100906A, such similarity of temporal 
behaviour is absent. It has a perfectly smooth optical light 
curve, gradually rising and decaying, a different profile from 
that in the 15-150 keV band (Fig.[3|). 

From this point of view, GRBs 100901A and 100906A 
support the earlier noticed diversification for optical emis- 
sion detected during the first few minutes after the onset 
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Table 14. Parameters of the time fits to the GRB 100906A Ught curves, where slope a is quoted in Fu oc t~°'v~^ . The first three fits 
are done with the broken power-laws, and the last one, with the power law. Also given are the spectral indices in the XRT spectral range 
(/3xrt)i retrieved from the Swift Burst Analyser and averaged over the relevant time intervals. The time intervals are the two parts of 
the whole time span (2nd column), divided at the break time (5th column). 



Light curve 


Time (ks) 


before break 

o (/3xrt) 


^brcak (ks) 


after break 
a (/3xRT> 


xVdof 


3-10 keV 
R (DSN) 
R (OSN+NOT) 
1-3 eV (MASTER) 


10 - 300 
28 - 54 
28 - 310 
0.18 - 2.52 


1.68 ±0.02 
1.36 ± 0.05 
0.14 ±0.02 
1.08 ± 0.02 


1.0 ±0.3 

1.1 ±0.2 
1.1 ±0.2 
1.3 ±0.4 


51+3 


2 q+°i 
wtf 

9 i7-t-0.03 
^•-^ ' -0.04 


1.1 ±0.2 

1.0 ±0.2 


145/78 
56.0/55 
135/58 
13.9/12 
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Figure 10. Spectrum of GRB 100906A compiled from the XRT 
data in the interval 45-65 ks and optical BVIR observations ob- 
tained by OSN around 14 h after the trigger (in the interval 50.9- 
54.0 ks) and corrected for the Galactic extinction. Solid line is the 
best-fitting absorbed broken power law model, whose parameters 
are listed as Fit 100906.4 in Table [TS] 




1000 



Figure 11. GRB 100906A light curves: unabsorbed 0.3-10 keV 
flux (black dots) observed by Swift/XRT, unfiltered MASTER 
band flux (blue crosses), i?-band flux by the OSN (red dots), and 
R-band flux by the NOT (green triangles). The optical points 
are corrected for the Galactic extinction. Broken-power law fit is 
plotted over the 0.3-10 keV light curve in white and black color, 
before and after the break time at 51 ks, respectively. Thin-line 
rectangles show the time intervals selected for spectral analysis. 



of a GRB (e.g., IVestrand etHlbOOSl . |2006l ): in some cases 
the prompt optical flux variations are correlated with the 
prompt 7-emission, and in other cases early optical after- 
glow Ught curve is uncorrelated with the 7-ray hght curve. 

Such division of GRBs can be understood within the 
picture of the optical emission generated from different sites 
of the relativistic jet: i n the internal, external forward, or re- 
verse shock waves (e.g.. lSari fc Piranll 19991 : iMeszaros fc ReesI 
I1999I ). The external forward shocks arise due to the interac- 
tion of the ejected material with the ISM and propagate into 
the ambient media, the reverse shock propagates into the 
ejecta itself, and the internal shocks are produced because 
of the collisions among many engine-eje cted shells (e.g., 
iMeszaros fc Rocs' 1 9971 : i Sari fc Puranll 19971 ). The correlation 
between the optical, 7, and X-rays can be achieved, if the 
optical emission is produced by the synchrotron mechanism 
within the same region as the hard radiation. 

The early X-ray and soft 7-ray spectrum of 
GRB 100901A shows a dependence close to which is 

typical for the synchrotron radiation ii < v < Vrn (where 
I'm and Vc are the characteristic synchrotron and cooling 
frequencies, respectively) in t he regime when ele ctrons lose 
energy very quickly (see, e.g.. lSari. Piran fc Nara van 19981; 
ISari fc PiranI Il999l ). Adding the MASTER optical point 
(corrected for the Galaxy extinction), and fitting simulta- 
neously the optical, 0.3 — 10 keV, and 15 — 150 keV data by 
an absorbed power law, we can determine the magnitude of 
the optical extinction in the host. Consistent values of 
and A^^' are obtained for the prompt spectral distributions 
(Table nU)) . indicating a persisting relation between the op- 
tical and high energy part of the spectrum. Basically, this 
means that the optical emission is the low-energetic pho- 
tons from the same region that generates the high energy 
photons. 

We should add that the data can also be reconciled with 
the prompt spectral distributions of GRB 100901A peaking 
between the optical and X-rays, if Vo-pt < I'c < v^- Then, 
a migrat ion of the cooling frequency Vc (oc t'^^"^ , see, for 
example. ISari. Piran fc Naravanlll998l ) may be hidden in the 
errors of resulted Ay*. Spectral fitting of an absorbed power 
law to the late-time MASTER and 0.3-10 keV data provides 
A'^* = 0.45±0.15, lower than A'^* = 0.8±0.1 for the prompt 
fits (values are for the SMC type of the extinction law; see 
i) l4.2p . This could be a confirmation of the above suggestion, 
but a detailed discussion is beyond the aims of the present 
study. 

In contrast to the above picture, GRB 100906A exhibits 
no correlation between the optical and high energy emission 
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at 15-150 keV. We suggest that in this case the optical ra- 
diation originates at the front of the external shock, and 
the light curve shows a smooth conversion of a bow shock 
to a self-similar mode of the afterglow with a power de- 



cay Jopt 



ex t 



with Q = 1.08 ± 0.01 (Table [H see also 



iKuvshinov et al.ll2010l ). 

A critical test for determining the site of the optical 
emission could be the observation of polarization. If radi- 
ation in the internal shock occurs in the order ed magnetic 
field, some polarization is expected (see, e.g., Isteele et al.l 
120091 '). This imphes that the magnetic field coherence scale 
is larger than the size of the visible emitting region. Thus, 
a discovery of a considerable (> 10%) polarization of the 
optical radiation, correlating with the prompt 7-ray emis- 
sion, would be a strong evidence in favor of the opti- 
cal radiation originating in the internal shock. Conversely, 
the afterglow radiation from the external shock occurs in 
the chaotic field of the interstellar medium compressed 
by the shock and should not have a considerable polar- 
ization iGruzinov &: WaxmanI [T999I ) . Low degree of polar- 
ization of the optical light would rule out the presence 
of a large-scale ordered m agnetic field in the emitting re- 
gion (|Mundell et al.ll2007l '). 

Fig. U shows our simultaneous observations of 
GRB 100906A with two orthogonal polaroids. The signals 
in the two channels are the same with an accuracy of 0.5 
per cent. Absence of the strong polarization is consistent 
with the picture of the different producing sites of the high- 
energy and optical radiation in GRB 100906A. However, we 
should note that the MASTER polarization observations of 
GRB 100906A were made with two optical tubes of one tele- 
scope. Such mode of observations is restricted in the sence 
of giving definite information on the linear polarization (see 
the Appendix). Polarization filters are oriented differently in 
the celestial coordinate system at the telescopes of the MAS- 
TER net. A future observation of a suitable GRB by two 
MASTER telescopes should thus provide an accurate polar- 
ization information. For the observations of GRB 100906A, 
we can infer that in the time interval from 100 s to ~ 1 h 
the degree of linear polarization of the optical emission is 
less than 10% with probability of about 60 per cent (see 
the Appendix and Fig. I14p . The full information of the in- 
ferred fractional linear polarization can be retrieved from 
Fig. [21 We note that the probabilities are obtained under 
a very conservative assumption that a source can have any 
fractional linear polarization from to 100%. 



5.2 General spectral properties 

In order to view spectral characteristics of the observed 
GRBs at one glance, we construct in Fig. [12] a diagram of 
the optic al flux density versus X -ray flux density at 3 keV, 
following I Jakobsson et al] l|2004 ). However, we plot not the 
values extrapolated to 11 h time mark, but rather the evo- 
lution for each GRB. Here, a spectral index /3 (Fi, oc i^"'') is 
shown by straight lines. Each optical observation is corrected 
for the Galactic extinction at a speciflc waveband using the 
NED Galactic Extinction calculator. As the intrinsic spec- 
trum is not known a-priori, we do not include a correction 
for the host extinction. The 3 keV flux densities are the best- 
fltting values obtained with an absorbed power law model in 



the XRT/ Swift energy band. For GRBs 100902A, 100905A, 
and 101020A, only upper optical limits are shown. 

The hmit for GRB 100902A at 132 - 162 s yields that 
the burst was rather X-ray bright and very faint in the op- 
tical waveband at the time, possibly due to a large intrin- 
sic ex tinction. Provided that T90 = 429 s (jStamatikos et ahl 
|2010| ). one can expect that the GRB would have low /3ox 
at (0.3 — 0.4) T qq. It might be a 'dark' GRB according to 
classiflcation by jjakobsson et "all ((20041 ) . bearing in mind its 
very prominent position on the diagram. GRB 100902A ap- 
p ears to be darker tha n any of the GRBs of the sample of 
Jakobsson et al.l l|2004l . see their Fig. 1), with all of them 
being above the /3ox ~ line. 

An anti-correlation between excess absorption column 
densities and redshifts of GRBs is expected due to the 
shift of the source rest-frame energy range below 1 keV 
out of the X - Ray T elescope observable energy range. Thus 
iGrupe et all (|2007l ) suggest a relation to estimate an up- 
per limit on the redshift. From the late-time spectrum of 
GRB 100902A, one finds an excess absorption column den- 
sity of ~ 2.8toi X 10^^ cm"2 and upper limit for z « 3.2 
(substituting lowe r value of A^h from expression (1) of 
|GruEeeEalJ(|2003)). High redshifts are thought to be pos- 
sibly responsible for darkness of GRBs, but the redshift ob- 
tained above is not h igh enough for this scenario (see, e.g., 
[Jakobsson et al.ll2004|. and ref erences therein); on the other 
hand, expression bv lGrupe et a l. (2007) is inferred for non- 
dark GRB s and might be u nappl icable. We also note that 
analysis of ICampana et all ((20101) who use a cutoff power- 
law model suggest a higher redshift, z ~ 4.5. 

The MASTER observation of GRB 100905A was done 
at the time interval 150-180 s that is much later than its 
Tgo ~ 3.4 s. I t can belong to the class of intermediate du - 
ration GRBs ((Horvathlll99i . [20021 : [Mukheriee erai1[l998l '). 
At 150-180 s after the trigger its XRT spectral slope is 
/3x = 0.2 ± 0.2. If the optical emission belong to the same 
spectral distribution, it would explain its weakness, because 
the line /3ox ~ 0.2 intersects the region of the 100905A on 
the diagram at the _Fopt values well below the optical limit. 

GRB101020A could not be observed by Swift/XRT. 
The fiux density at 3 keV for time 0.65 x Tgo, where T90 = 
177 s, is extrapolated from the BAT data and has a consid- 
erable uncertainty, thus resulting in a wide area in Fig. 1121 
Nevertheless, the optical-to- X-ray flux ratio can be similar 
to that of GRBs 100901A and 100906A. 

We see from Fig. [T2| that the flux density ratios of 
GRB100901A and GRB 100906A ultimately tend to the 
area between (3 = .5 and 1.25. This is consisten t with the 
fireball model (e.g. [Sari. Piran fc Naravarj [l998h . Spectral 
slope ~ — 1 at late times is consistent with the prediction 
for the slow-cooling regime for frequencies higher than the 
cooling frequency for electrons with distribution of Lorentz 
factors N{-ye) oc 7"^ with p > 2. 

The distinct prompt behavior of GRB 100901A and 
GRB 100906A is evident in the diagram. We witness the 
optical and X-ray fiux density correlation in GRB 100901A 
at the time of the last 7-ray pulse, together with an ap- 
proximately constant spectral index. In GRB 100906A, dur- 
ing the 7-ray pulse, occurring about 100 s after the trigger, 
the 3 keV flux density undergoes a strong variation. The 
same can be said about the spectral indices in the XRT and 
BAT bands. However, the optical flux density is roughly at 
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Figure 12. Diagram of the optical flux density versus X-ray flux density at 3 keV for five GRBs observed by MASTER in September 
and October 2010. Each point represents a simultaneous observation in the optical and by 5toi/t/XRT. The optical flux densities have 
been corrected for the Galactic extinction. Diamonds connected with grey line show evolution of ORB 100901A (eight time intervals: 
193-233, 312-372, 387-467, 481-581, 595-715, 3800-4150, 27000-28500, and 102800-121500 s); squares connected by blue lines, of 
GRB100906A (seven time intervals: 84-94 (point 1), 105-125 (point 2), 330-390, 410-480, 942-1112, 6000-7000, and 10300-11000 s). 
For GRBs 100902A, 100905A, and 101020A, only upper limits on optical flux density are given at post-trigger times 132-162, 150-180, 
and 106-126 s, respectively. When tracks are parallel to /3 =const lines, the optical flux density is emphatically correlated with the 
3 koV flux density We adopt Galactic extinctions Ay = 1.037 for GRB 100902A, Av = 0.191 for GRB 100905A, and Av = 0.051 for 
GRB 101020A from the NED Extragalactic Calculator. 



the same level. We remind that the spectral fitting provides 
an evidence that the prompt high energy SED has a peak, 
which migrates from about 30 to 4 keV at times 80-130 s 
after TO (see § 14. ip . It should be added that during those 
times the X-ray emission in the interval of 0.3-10 keV can 
be contributed by the two different sites with a proportion 
varying with time: the one that is responsible for the 7-rays 
and one producing the optical photons. 



5.3 Amati and Ghirlanda relations for 
GRB 100901A 

The total duration BAT spectrum of GRB 100901A yields 
Eiso = 6.3 X lO''^ erg in the energy interval 1 - 10 keV 
in the rest frame. From the Amati relation (|Amatil [2OO6I ') 
one gets the expected position of of the energy distribution 
peak in the rest frame Ep^i — 230 keV, which translated to 
the observer frame becomes Ep ~ 100 keV. This does not 
contradict with the lower limit i?™'" ^ 75 keV, derived on 
the basis that no characteristic energy can be found form the 
spectral analysis of the BAT data, and the photon index of 
the power law is about -1.5. 

The 0.3-10 keV light curve oc t"^ '^ from 4 x 10* s to 
about 10^ s when a possible steepening occurs. The MAS- 
TER W-band hght curve enables only a rough evaluation 
of temporal index in the interval 3 x 10'* — 3 x 10^ s: it is 



decaying approximately as ~ t~^'^ with just one late op- 
tical point at ~ 10^ s. Using the tentative jet break time 
~ 10® s, one gets the opening half-angle of the jet ~ 13 deg 
for efficiency of converting the ejecta energy into 7-rays 
rj-y — 0.2 and the number density of the ambient medium n — 
3 cm ^. Few jets are reported to have comp arably wide an- 
gles (|Ghirlanda. Ghisellini fc Lazzatill2004l ). Corresponding 
coUimation-corrected energy would be Ej = 1.7 x 10^^ erg, 
which using Ghirlanda relation predicts _Ep,i ~ 700 keV, a 
value about 3 times higher than prediction of the Amati 
relation. 

We conclude that GRB 100901 can be consistent with 
the Amati relation, if Ep^i ~ 230 keV, but the Ghirlanda 
relation is hardly applicable. 



5.4 Amati and Ghirlanda relations for 
GRB 100906A 

iGolenetskii et ahl (|2010l ). fitting the spectrum of 
GRB100906A in the time interval 3 - 142 s with the 
Band function, derive values -Epcak ~ 142lg|!|^ keV 
and _Biso = (2.2 ± 0.4) x lO'^^ erg. The peak of 
the energy distribution in the rest frame is at 
£p,i = {1 + z) X Epcak = 387l?g4 keV. Central val- 
ues of Eiso and Ep^i are cl ose to the A mati correlation 
£p,i = 95 X Ef^f (figure 2 of lAmatillioO^ . 
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If we assume jet break time % = 51 ± 4 ks (see 
§ 14.41) and substitute the above -Ejso into equation (1) 
of lGhirlanda. Ghisellini fc Lazzatl (|2004l l (where we assume 
the values of the efficiency of converting the ejecta en- 
ergy into 7-rays -q-y = 0.2 and the number density of 
the ambient medium n = 3 cm^), we obtain the open- 
ing half angle of the jet = 3.31 ± 0.08 deg. Then 
the coUimation-corrected energy is Ej ~ (3.7 ± 0.7) x 
10^" erg. The last value together with Ep^i agree with the 
Ghirlanda relation gp.i ~ 480(-E-y/10 '^^ erg)'' '^ (figure 1 of 
iGhirlanda. Ghisellini fc Lazzatill2004l ). This is especially so 
due to the big uncertainty in the observed -Epeak value of 
GRB 100906A. 



5.5 Long life of a central engine 

A prolonged activity of the central e ngine was suggested a s 
a source of an extended afterglow bv lKatz fc Pirai] fTooJ). 
It was considered later as an explanati on for plateaus on 



X-ray light curves o f some GRBs (e.g. 



ray 

200ll: IZhang et al.l l2006l: 



Zhang fc Meszaroj 
Liang. Z hang fc Zhand l2007l : 



Willingale et al.ll2007l : TLipunov fc Gorbovskov 2008 ) . In par 



ticular, sharp drops at late times on the light curves indicate 
that suc h variations are l ikely due to some internal mech- 
anisms (|Troia et al.l l2007l : fzhang 200i). To account for a 
lasting activity of the central engine of GRBs, a model of a 
rapidly-spinning, highly-magnetized neutron star, a 'proto- 
magnetar', formed in the result of a core collapse, was pro- 
posed (see review by iMetzgeJ I2OI0I '). Another object that 
can provide a long-term energy supply is a rapidly-spinning, 
magnetized collapsed core with mass exceeding an upper 
limit for a neutron star. Its fast rotation stops the collapse, 
and a quasi-stationary configuration is pre sumably formed 
at this stage, which is calle d a 'spinar' (|Lipunoval 1 19971 : 
iLipunov fc GorbovskovlboOTi '). Before a collapse to a black 
hole is possible, enough angular momentum should be taken 
out from the system. The time required for angular mo- 
mentum dissipation determines the duration of the central 
engine activity and depends on the magnetic field strength, 
as the energy and angular momentum losses are governed 
by t he magneto-dipole mechanism . 

iLipunov fc Gorbovskovl (|2008l ') proposed a set of 1-D 
equations determining evolution of a collapsing rotating 
magnetized core with main relativistic effects taken into ac- 
count. These equations describe the principal behavior of 
both the proto-magnetar and the spinar. Solutions suggest 
two peaks in energy release: the first takes place when dy- 
namical collapse is halted by centrifugal or pressure forces, 
and the second, when final contraction occurs. The peaks are 
separat ed by a plateau. Three exempl ary solutions of equa- 
tions bv lLipunov fc Gorbovskovl (|2008l ') are shown in Fig. ll3l 
The first peak is set to occur at zero time and not shown in 
detail. In the model, the first peak is evaluated on a dynami- 
cal time-scale and does not involve specific jet and afterglow 
radiation mechanisms. If the first stage of a collapse finishes 
close to Tg — G M/(? , then the second peak is extinguished. 
A power of a proto-magnetar, defined by the pulsar dipole 
mechanism, approaches the t"^ law after the second peak. 
A power of a spinar has always a very sharp drop after the 
second peak. For each of the models in Fig. 1131 calculated 
power is multiplied by an arbitrary factor of 0.1. 

Two curves with lower plateaus represent models with 
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Figur e 13. Energy release o f a col lapsing core calculated fol- 
lowing ILipunov fc Gorbovskovl l|2008h . Points with bars arc the 
obs erved 0.3-10 keV and 15-150 keV flux f rom GRB 100901A 
(Sa kamoto et al.ll2010bl : IPage fc Immlei]l201oh translated to red- 
shift 1.408. Three curves with plateaus represent numerical solu- 
tions for energy release rates, multiplied by 0.1, for cores of dif- 
ferent mass and angular momentum. The dashed and dot-dashed 
lines are solutions for proto-magnetars of the same mass and with 
different angular momentum at the halt of the initial dynamical 
collapse that produced the first peak before the plateau. The first 
peak is referred to time zero; the details of the first peak are of no 
relevance. The top solid line shows a solution for a collapsing core 
with such a mass that makes a collapse to a black hole inevitable. 
For specific parameters of the models, see Discussion. 



mass M = 1.9 Mq. They both go as oc t~ at late times. 
In the model with monotonic transition to a decay, a proto- 
magnetar forms with period P = 2 ms and poloidal mag- 
netic field B ~ 4 X 10^* G at ~ 13 rg. The middle plateau 
represents a newly born magnetar of higher angular mo- 
mentum: the dynamical collapse halts at about 40 rg with 
P = 5 ms and B ~ 7 x 10^^ G. The object contracts fur- 
ther, to 10 Tg, causing a spin up to P = 1.7 ms and the 
dipole magnetic field amplification up to ~ 10^^ G. This 
specific evolution manifests itself as a bump at the end of 
the plateau. The solid curve shows a solution for a spinar 
with M = 4 Mq . At the plateau, the dimensionless Kerr pa- 
rameter OKerr = Jc/GM^ = 3, P ~ 1 ms, B ~ 1 .5 x 10" G, 
and the characteristic radius 15 rg. The spike at the end is 
numerically solved and indicates the collapse to a black hole. 
The initial magnetic field strengths are chosen in order to 
obtain the characteristic time of a plateau ^ 10'* s. 

A spinar and a proto-magnetar, both can explain long 
activity of GRBs. The law of temporal decay after the second 
peak (or the end of a plateau) can be defined by an afterglow 
and does not allow us to decide which compact object is 
behind the scene. The efficiency factor, which is set 0.1 for 
Fig |13l cannot be constrained in our general approach. 



6 SUMMARY 

The net of optical robotic telescopes MASTER coming into 
full operation over the Russia territory is encouraged by the 
need of dense early observations of GRBs around the globe. 
Optical ground-based observations, with polarization mea- 
surements, are indispensable in revealing the nature of the 
GRB central engine. The MASTER team is working on ob- 
taining optical linear polarization information at the earliest 
time since a GRB trigger, before 7-emission ceases, by ob- 
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serving simultaneously with several telescopes. Our goal is 
to deliver full polarization information for GRBs in the op- 
tical band as early as possible after the trigger time. 

Here we presented the early results of the MASTER net- 
work. During 2010, MASTER observed five GRBs described 
in this paper. Optical transients were found for two bursts: 
GRBs 100901A and GRB 100906A. Their detailed analysis 
was done basing on the MASTER data, the data from the 
1.5-m telescope at Sierra Nevada Observatory and 2.56-m 
Nordic Optical Telescope, as well as the publicly available 
data from the Swift mission. 

The light curves and spectra of GRBs 100901A and 
100906A suggest different origin of the prompt optical emis- 
sion. The correlation between the optical and high-energy 
light curves of GRB 100901A, as well as the results of spec- 
tral fitting favor the hypothesis of the common origin for the 
optical and high-energy spectrum. During the 7-ray flare 
started at about 300 s after the trigger, the slope of the 
spectral distribution from the optical to 150 keV is close to 
— 1/2. This is expected for the synchrotron radiation at the 
frequencies above the cooling and below the characteristic 
synchrotron freque ncy in the 'fast cooling' regime (see, e.g. 
loan fc PiranllT999ll . The resulted values of the host optical 
extinction are consistent for the prompt fits and similar to 
the values obtained by fitting an absorbed power law to the 
late-time spectral data. This means that the prompt optical 
points lie on the same spectral distribution which describes 
the 7- and X-ray data, supporting the hypothesis of their 
common production site and mechanism. The estimated ra- 
tio of the GRB host hydrogen column density A*'^* to the 
total extinction J^y^ suggests that the host of GRB 100901A 
possesses the metal-to-dust ratio comparable to the level 
measured for the Milky Way or Magellanic Clouds. 

In GRB 100906A, no apparent correlation between the 
optical and 15 — 150 keV light curves is seen. We suggest that 
in GRB 100906A the prompt optical radiation originates in 
a region different from the 7-ray production site, probably, 
in a front shock. An analysis of GRB 100906A reveals that 
a spectral component, responsible for the 7-ray pulse at ~ 
100 s after the trigger and the downward migration of the 
X-ray peak on the spectral energy distribution, does not 
correlate with the optical emission at these times. 

Signals in two orthogonal polarizations, measured by 
the MASTER telescope at Tunka from GRB 100906A, are 
found to be equal within 0.5 per cent. We calculate the prob- 
ability distribution for the fractional linear optical polariza- 
tion of GRB 100906A. We point out that a future observa- 
tion of a GRB by at least two MASTER telescopes should 
yield an accurate polarization information. 

The 0.3-10 keV light curve of GRB 100906A has a break 
at around 14 hr after the trigger. The X-ray spectral index 
(non-variable around the break) and the temporal slopes 
can be reconciled with the closure relations for the rela- 
tivistic jet in the 'slow cooling' regime before the break 
and the spreading uniform jet after the break. It is not 
clear from the data in the R band if this break is achro- 
matic. Estimated jet opening angle, 3.31 ±0.08 degrees, and 
coUimation-corrected energy, E-i ~ (3.7 ± 0.7) x 10^" erg, 
imply that GRB 100906A follows the Amati relation and 
does not contradict the Ghirlanda relation. An analysis of 
the late-time spectra provides the upper limit on the optical 
extinction in the GRB host. 
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Appendix. Estimating degree of linear polarization 
from observations with two fixed polarizing filters 

Observations of bright GRB 100906 have been made with 
just one telescope able to provide polarization measurements 
(the MASTER telescope at Tunka). Generally, if an obser- 
vation is made only with two polarizing filters at a single set 
of positions, one can obtain just an estimate on the degree 
of linear polarization. 

In the consideration below we assume zero circular 
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Figure 14. Left axis: likelihood C{Pl > p™'" | ^ /j), where 
Ii,2 are the signals from two orthogonal polarizing filters, is 
an unknown fractional linear polarization of a source, and p™'" 
is a value plotted along the x-axis. The several likelihood func- 
tions are plotted as solid curves marked with the values of o-p, 
the known relative error between the measurements with two or- 
thogonal polarizing filters. The dashed line shows the probability 
of an unobserved source to have Pl ^ pmm mjjgj. ^j^g prior hy- 
pothesis that a source has no preference to any degree of linear 
polarization. Right axis: the ratio of the total intensity / to the 
measured signal Ii for \D\ < cTp = 2%. Its limits are shown by 
the double dotted line. 

polarization. Let Pl be the fraction of linearly polar- 
ized radiation and 6 is the angle of maximum polariza- 
tion, i.e. the angle between the a;-axis of the reference 
system and a segment which a polarization ellipse de- 
generates into for the case with zero circular polarizatio n 
(see, e.g.. lLandi Degl'Innocenti. Bagnulo fc Fossatl l)2007| )). 
Then Pl = ^/Pg + Pd, where Pq = Q/I, Pu = U/I, and 
I, Q, U axe the first three Stokes parameters; / is the to- 
tal intensity with no polarization filter interposed (in this 
section we do not consider color filters). 

To derive value of Pl one needs to perform observa- 
tions with filters for linear polarization positioned at three 
angles, for example, at 0°, 45°, and 90° with respect to the 
reference direction. We have two orthogonal filters at each of 
the MASTER-II telescopes; thus, if such telescope observes 
alone, only a lower limit on Pl can be deduced. If this lower 
limit is a virtual zero, there is still some probability that 
actual linear polarization is substantial. 

Let Ii and I2 be the signals of the detector with filters at 
0° and 90° interposed. Then a value D = {h- l2)/{h + h) 
can be calculated, which is expressed as follows: 

D = Pl cos 26 . (1) 

If the value of D differs from zero {\D\ is greater than 
the relative error Up between measurements with two dif- 
ferently positioned filters), a certain amount of linear polar- 
ization is present in the hght with a lower estimate for it: 
pmm = |/)|^ calculated as if a filter is interposed at angle 0° 
or 90° with respect to the direction of polarization plane of 
the incident light. 

If the module of D is less than o-p, then formally we 
have zero as a minimum estimate for the degree of linear 
polarization. 100-percent linear polarization is also possi- 
ble in this case, if a filter is accidentally positioned so that 
a source's polarization plane is inclined by 45° to it. This 



possibility can be numerically evaluated considering a to- 
tal range of the filter positions that would result in a value 
of |D| < (Tp. Substituting D ^ ap and 9 = 7r/4 + S8/2 in 
eq. IT]), for any preset value of Pl we get a value of the 
allowable interval 59. As our configuration consists of two 
orthogonal filters, source polarization plane should be in 
either of four intervals 66 around angles A'^ x 7r/4, where 
A^ = 1, 2, 3, 4. Thus, a source with linear polarization degree 
Pl yields value \D\ < ap with probability 466 /2n. Now we 
can calculate the conditional probability of a source with 
observed \D\ < ap to have a particular degree of linear po- 
larization Pl using the Bayes' theorem. We make a conserva- 
tive assumption that a source can have any fractional linear 
polarization from to 100%. From this, the likelihood func- 

tion c{Pl > pr" ! Ji ^ 12) = £.{Pl > pr" I 1^1 < ^p) 

can be derived, where Ii f» /2 is a realized event, Pl is the 
actual fractional linear polarization of a source, and P™'" is 
some value of it. The probability that the actual linear po- 
larization of a source exceeds value Pl is shown in Fig. 1141 
Thus, for example, if no difference between Ji and I2 is reg- 
istered in the measurements with 2 per cent error, there is 
a fifty-percent chance that a source has linear degree of po- 
larization greater than « 6 per cent. 

Accordingly, a total intensity from the source is greater 
than observed 7i ~ l2- One can see from Fig. [14] that in 
the observations with |D| < ap the total intensity is always 
about the same factor greater (twice) than measured 7i. This 
can be seen considering expression for the signal from one 
polarizer: h ^ I {1 - Pl)/2 + I Pl cos^ 6. For small Pl the 
light behaves as unpolarized, and for large Pl the allowable 
interval S6 approaches zero and 6 ~ 7r/4. 



